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Abstract 

We investigate the continuum limit in SO{N) Principal Chiral Models concen- 
trating in detail on the 50(4) model and its covering group SU{2) ® SU{2). We 
compute the mass gap in terms of A-j^ and compare with the prediction of Hol- 
lowood [1] of m/ A-j^ = 3.8716. We use the finite-size scaling method of Liischer 
et al. [2] to deduce m/ A-j^ and find that for the 50(4) model the computed result 
of m/ A-j^ ~ 14 is in strong disagreement with theory but that a similar analysis 
of the SU{2) (8) SU{2) yields excellent agreement with theory. We conjecture that 
for 50(4) violations of the finite-size scaling assumption are severe for the values 
of the correlation length, ^, investigated and that our attempts to extrapolate the 
results to zero lattice spacing, although plausible, are erroneous. Conversely, the 
^ finite-size scaling violations in the SU (2) Cg) SU (2) simulation are consistent with 

Q perturbation theory and the computed /3— function agrees well with the 3-loop ap- 

proximation [2] for couplings evaluated at scales L/a < ^, where ^ is measured in 
units of the lattice spacing, a . We conjecture that lattice vortex artifacts in the 
50(4) model are responsible for delaying the onset of the continuum limit until 
much larger correlation lengths are achieved notwithstanding the apparent onset 
JnI of scaling. Results for the mass spectrum for SO(N), # = 8, 10 are given whose 

comparison with theory gives plausible support to our ideas. 
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manifold M = G/H^ where G and H are Lie groups and H C G, has provided a 
deep understanding of field-theoretic methods applied to both high energy physics and 
critical phenomena. In particular, the continuum field theories derived from these 
models are renormalizable and asymptotically free and so act as good models to test 
ideas applicable to four-dimensional gauge theories. For many models the exact S- 
matrix for the continuum field theory has been conjectured [3, 4, 5] and this allowed 
a number of workers [5, 6, 7] to give exact predictions for the mass spectrum and 
to use the Thermodynamic Bethe Ansatz together with the property of asymptotic 
freedom to predict the ratio ni/Aj^ for these theories, where m is mass-gap and 
sets the mass scale in perturbation theory. It is important to note that solely the 
local properties of M, i.e., the algebra which generates infinitesimal displacements on 
M, and not the global properties of M are used to obtain these results. To be able 
to control the approach to the continuum limit for a lattice field theory is central 
to current calculations in lattice QCD where criteria based on scaling are used to 
determine whether dimensionless ratios for computed observables have attained their 
continuum values and so can be compared with experiment. The existence of theoretical 
predictions for 2D chiral field theories which link the short and long range structure of 
the theories gives a unique opportunity to study the approach to the continuum in a 
lattice model and also to test the assumptions underlying the theoretical analysis. 

The correlation length is most conveniently measured in units of the lattice spacing, 
a, and is denoted by the dimensionless variable, ^, and the mass-gap is defined by 
m = For ^ sufficiently large the mass spectrum inferred from the correlation 

functions of suitable chosen interpolating operators and the mass-gap measured in units 
of Ajyj^ should agree with theory. Such agreement has been found for the SU{N) matrix 
models [8, 9] for values of = 3 — 15 already for moderately small correlation lengths: 
^ > 5. However, the manifold of SU{N) is simply connected but this is not so for many 
other models for which it is important to see if any obstructions/lattice artifacts delay 
the onset of the continuum limit. We may then learn how to be more confident in 
controlling the approach to the continuum. Similarly good agreement with theory has 
also been found by Wolff [10] for the 0(4) and 0(8) spin models but not for the 0(3) 
spin model [11, 12]. A comparison of 0(3) and RP^ models has recently been the 
subject of much attention [13, 14, 15, 16, 17] since they are expected to correspond to 
the same continuum theory although they differ in the global properties of the defining 
manifolds. It has, however, proved impossible to demonstrate the equivalence of these 
two models in the continuum by computer simulation alone. The study reported in this 
paper highlights the difficulty of reliably controlling the approach to the continuum in 
a number of common lattice models. 

We study the SO{N) matrix models on the lattice with the action given by 

S{V) = (3Lj2Tr{u{x,t)V^{x,t+l) + V{x,t)V^{x + l,t)) , (1.1) 

where U(a;, i) is an A X A orthogonal matrix. Our analysis mainly uses the finite 
size scaling method described by Liischer et al. in reference [2] which is based on the 
coupling introduced by Nightingale [18]. For the lattice regularized theory we define 
the coupling constant u by 

u{L/a,h) = m{L/a,h)L , (1.2) 
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length is determined by the decay of the correlation function on the L X T strip in the 
time direction. The physical properties can be separated from the cutoff dependence 
and the continuum coupling, u{Lmoo), can be defined by writing 

u{L/a,(3L) = u{Lmoo,a/L) = u{Lmoo) + 0{{a/LY) . (1.3) 

The physical scale is set by moo = l/f^^ooj the mass-gap in the infinite- volume limit. 
Apart from a trivial reparametrization, this decomposition relies on the physical hy- 
pothesis of scaling with the exponent u parametrizing the corrections to scaling. For 
sufficiently large /3/, and sufficiently small L/a perturbation theory can be used to 
compute M, and since it is proportional to 1//3l at tree-level it can be viewed as a 
definition for a running coupling of the continuum field theory. In particular u will 
evolve according to the universal part of the /3-function. 

u{Lmoo) can be calculated from a simulation of the lattice theory. For Lrrioo > 
1 this is straightforward: for a given value of /3/, compute arrioo and then simulate 
on a lattice of width L given in units of a by L/a = Lnioo/amoo, and so compute 
u{L/a, (3l) from equation (1.2). Repeat for a sequence of (3l values, the largest of 
which is determined by computer resources. Check whether the results are compatible 
with the scaling hypothesis (equation (1.3)) and, if so, take u{Lmoo, a/ L) with the 
largest L/a as the best approximation to u{Lmoo) or, better still, use equation (1.3) to 
fit the a-dependence and extrapolate to a = 0. 

However, this direct approach cannot be used for small values of Lrrioo since L/a > 1 
and arrioo is bounded from below by computer resources and so limits how small Lrrioo 
can be. This problem is resolved using the renormalization group by computing the 
lattice step-scaling function S: 

S(s,M(L/a,/3L)) = u{sL/a,(3L). (1-4) 

Since no reference is made to moo we can choose, for any given value of m, L/a as large 
as the computer resources will allow by tuning /3/, appropriately. Again we expect the 
scaling hypothesis to hold: 

S(s,M(L/a,/3L)) = (j{s,u) + 0{{a/LY). (1.5) 

(t(s, u) is the continuum step-scaling function determined implicitly from the /3-function 

by 

The /3— function has the perturbative expansion 
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[i{u) = -u'J2ku' , (1.7) 
/=o 

Having computed u{2), say, by the direct method a sequence of values m„ = m(2s") can 
be computed using 

Un = a{s,Un-l) . (1.8) 

A reasonable choice for s is s = 1/2. For couplings, m„, evaluated at sufficiently 
large n a{s,u) is well approximated by perturbation theory and A^ is given to 2-loop 
order by 

LA= (5o^*(Lmoo))-'^/'"expf--— ^ -) , (1.9) 

V hQu{Lm oo)) J 
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to overall rescaling) and for some cases 62 is known [2] yielding a 3-loop estimate for 
A^. A 1-loop calculation then relates A^ to ^j^- 

Alternatively, the approach used in references [12, 10, 8, 9] can be employed where 
a bare lattice coupling is defined to be a running coupling evaluated at the scale given 
by the lattice spacing, a. Two candidates are 

UL{amoo) = — : the 'lattice scheme' , 

PL 

(1.10) 

UEiamoo) = • 'energy scheme' , (l-H) 

where C = nf, nj/N for spin and matrix models respectively, and raj is the number 
of degrees of freedom in the field. The coupling ue was first defined by Parisi in 
reference [19]. For ^00 ^ 1 both of these couplings should scale with a according to 
the perturbative /3— function thus allowing aA/, and aA^ to be computed from ul and 
Ue respectively according to equation (1.9) with a replacing L. Since ^00 is known this 
allows m/ to be estimated once the 1-loop calculations relating to A/, and 

A^ have been done. 

The main motivation for our study was the recent work by Hollowood [1, 20] in 
which he has given theoretical predictions for m/ for all > 3, so allowing us to 
compare our lattice calculations with theory. The comparison between the studies of 
the lattice and continuum models sheds light on the approach to the continuum limit 
in the lattice formulation and can provide verification of the assumptions underlying 
the theoretical analysis of the continuum theory. We calculate u{Lmoo, a/L) defined in 
equation (1.3) by simulation for various values of anioo and discuss the extrapolation of 
the results to a = and we repeat the analysis for the covering group of 5*0(4), namely 
SU{2) (SiSU{2). We also compare the results with an analysis based on the lattice- and 
energy-schemes for the couplings defined in equation (1.11). For the 5*0(4) model we 
find no evidence that finite-size scaling holds for the accessible range of couplings, but 
that the results can deceptively suggest that it does hold since the violations of scaling 
diminish only slowly with increasing ^ so giving the semblance that they are negligible. 
In contrast, for SU{2) ® SU{2) finite-size scaling holds with O(a^) violations and the 
value deduced for ni/Aj^ agrees well with the prediction. 

We also compare the computed mass spectrums in the SO{N) model for = 8, 10 
with theory. The interpolating fields for the vector bound states of SO{N) can be 
easily constructed and A^ = 8 is the smallest value for which there exists more than 
one such state. The results show a persistent deviation from theory for the range of ^ 
investigated. 

In section 2 we present the relevant 1-loop calculations which relate the various 
A— parameters; in section 3 we describe the simulation and measurement techniques 
including a method for variance reduction; in section 4 we give the simulation results 
and analysis for the study of the 5*0(4) matrix model; in section 5 we give the simulation 
results for the SU{2) ® SU{2) model; in section 6 we discuss the A^ = 6, 8 models and 
compare the mass spectrum and Aj^g, computed in the lattice and energy schemes, 
with theory; in section 7 we discuss our interpretation of all the results and how further 
studies may elucidate our findings. We also draw our conclusions. 
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The energy-scheme can be related to the lattice-scheme using the 1-loop calculation for 
{E) expanded in 1//3l = ^L- This calculation is familiar and we just give the result for 
general SO{N): 

{E) = {l-l^Tr{\J{x,t)\J^{x + l,t))) . 

N -I f N , , , 

1 + ^:777^ + ■■■] ■ 2.1.1 



8/3l V 32/3l 
The energy-scheme coupling, /3^, is defined as 



N - I 



8{E) ' 

and the ratio between A-parameters is then 

Ae (A 

where = (3l — A . 

A calculation similar to that described in reference [21] gives the result 

is related to Al by a 1-loop background field calculation in which the 2D theory 
on a strip is converted to a local ID theory which can be directly solved for the mass 
gap since there is an equivalent Schrodinger equation. A similar approach was used in 
reference [22]. To analyse the ID system it is sufficient to study the action 

S{W) = /3QMj2Tr {WtWt+i) . (2.1.5) 

t 

We represent W by 

W = exp(^-T) , 

where the generators Ti are defined so that 



Tr{T,T,) = . (2.1.6) 



Then for small cf) we have 



S{cl>) ~ + ... . 

t 

This defines the quantum mechanics of a particle constrained to the SO{N) group 
manifold with mass = Pqm/'^- The corresponding Schrodinger equation is then 

-l-B'Yi{W) = EiYi{W) . (2.1.7) 
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The Yi are group harmonics with / a generic label, and satisfy [23] 

B^Yi{W)+CiYi{yV) = . 

Using lo cx 1 and yi(W) cx Tr{W) we find Co = and C\ = Cp, the Casimir of the 
fundamental representation deduced from the generators. From equation (2.1.6) we 
find 
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m{L) 

N -^ ^ ^ 

4/3 • 

^Pqm 

We now calculate [iqM in terms of [ip to 1-loop order. The field variables U on the 
2D strip of width L and length T can be expressed in terms of fiuctuations about a 
background field which is constant across the strip and slowly varying in the 
coordinate along the strip. We write 

\]{x,t) = wy%5^(^'*) W^^ , (2.1.9) 

where (j){x^t) = cf){x,t) ■ T is the fiuctuation about W^. The definition of Wt is given 
by the relation 

X 

with Sf a positive symmetric matrix. Using equation (2.1.9) this corresponds at lowest 
order to the constraint 

Y,Hx,t) = . (2.1.10) 

X 

Because is slowly varying we can expand the effective ID action in terms of 
where = • T and 

To 1-loop order this equation can be implemented by choosing 

Wt = 1 and Wt+i = e"^' . (2.1.11) 

Substituting for U(a;, t) and from equations (2.1.9) and (2.1.11) into the SO{N) 
action, equation (1.1) we find 



X,t 



x,t 

Expanding to quadratic order, choosing = 2/ (3^ and using equation (2.1.6) we find 

S = Sq -\- Si , 
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Si = L/3i^-iAt • At 

t 



(2.1.13) 

The abbreviations cj) for (j){x^t) and for + 1) have been used and V is the 2D 

lattice Laplacian. It should be noted that there is no term linear in A^. This is due 
to the specific form of the decomposition in equation (2.1.9) and the identity 

Tr {{T,T, + T,T,}Tk) = . 

The absence of such terms simplifies the calculation since all contributions to the ID 
effective action, S'g||-(A), are simply given by 

5,ff(A) = {Si) . (2.1.14) 

Other parametrizations require an evaluation of the quadratic terms in (Sj). The 
average is with respect to the fiuctuation measure exp(S'o) taking into account the 
constraint in equation (2.1.10) which eliminates the zero mode in the s-direction. We 
use the Gaussian results 

{Ml) = ^(0,0)5,, {M'j) = G{0,1)S,,. (2.1.15) 

^(O, 0) and G{0, 1) are given by the expressions 

" 2^ ^1 4(sin2(7rp/L)+sin2(g/2)) ' 

^(0,1) = G'(0,0)-i. (2.1.16) 
We use generators defined in equation (2.1.6) with k = 1 and the identities 



N - 1 

^^^^ - 7 Tr{T^TjTk) = -^ftjk ■ 



(2.1.17) 



Then 



{Tr{<l>'A')) = {Tr{<l>"A')) = ^ G(0, 0) (-iA • A) , 
{Tr A')) = ^G(0,l)(-iA.A) , 
(rr(<^A<^'A)) = ^G'(0,l)(-iA.A) . (2.1.18) 
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Using equations (2.1.13) to (2.1.18), we find 

5eff(A) = L/3,ff^(-iA.A) , (2.1.20) 

t 

wliere 

/3eff = &-^(^logi--4) . (2.1.21) 

Note tliat tlie coefficient of log L is 5o for the coupling ul = 1//3l- Substituting into 
equation (2.1.8) with [iqM = -^/3gff gives 

1 _ 1 ^ 4/3,ff 
C TV - 1 

(/3l - 5o(log L + 7^)) , 



m(L) m(L)L 7V-1 

4 



TV - 1 

7jv = ^ 7 - 27rA . 2.1.22 

' 2 TV -2 ^ ' 



We thus deduce that 



exp (7 



N 



2{N 

Using equation (2.1.4) we finally get 

A- 



exp -27rA+— — . (2.1.23) 



''MS 

Ac 



\/32exp(27rA) , (2.1.24) 



which is independent of N . 

From equation (2.1.22) we also deduce the tree-level relation 



N -1 , , 

——UL . (2.1.25) 



2.2 0(iV)-spin models 



The results for the mass-gap for 0(A^)-spin models is given to 3-loop order in [2] 
and so we include only a brief outline of the 1-loop effective ID calculation here for 
completeness. A similar calculation was done by Liischer [24] but it is instructive to 
present it in a concise formulation consistent with the previous section. To relate the 
energy- and lattice-schemes we use the expansion [10] 

N -1 f 1 

and the definition 

/I - 
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From ref. [21] we have 

^ = V32expf , , 
Al ^V2(7V-2) 

To calculate m{L) on the L X T strip we use the same method as in the previous 
subsection with the 0{N) spins s{x,t) expressed in terms of the background field 
and the fiuctuation field cf){x,t) as 

s{x,t) = St^l - cf){x,t) ■cf){x,t) + cf){x,t) , (2.2.1) 
where • = 1 and cf){x, t) is an {N — l)-dimensional vector satisfying the constraint 

J2cf>{x,t) = . (2.2.2) 

This constraint eliminates the zero-mode divergence in the calculation. We then choose 

St = So and St+i = e^'^^^ Sq , (2.2.3) 

where Sq = (1, 0, 0, . . . , 0) and is the slowly varying ID background field. T12 is 
the generator of SO{N) rotations in the 12-plane. As in the SO{N) case the particular 
form of the parametrisation in equation (2.2.1) ensures that the ID effective action is 
given by 

•^eff = i^l) ' 

there being no linear terms in in the expansion of the S{s) in A and cf). The 
calculation follows the same steps as in the SO{N) case and gives 

/3eff = h - {N - 2) Q-log L - - \ , {2.2 A) 

where the coefficient of log L is identified with bo and, as before, A = 0.0351637. Then 



Using equation (2.2) we conclude that 

Ac 



y32exp(27rA) , (2.2.6) 



which is identical to the SO{N) result. The two results in equations (2.1.24) and (2.2.6) 
must be the same for = 3,4 since the two models have the same continuum limit. 
That they are independent of N ensures the results are identical for all N . 
The tree-level relation between couplings is then 

N -1 , , 

u = — —UL . (2.2.7) 
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ref. [25] we used an overrelaxation update [26, 27] applied to embedded 0(2) models. 
In terms of CPU-time requirements this algorithm out-performed a multigrid algo- 
rithm [28] for correlation lengths up to about 20 in the case of the CP'^-model in two 
dimensions. 

In order to save random-numbers, and hence CPU-time, a large fraction of the 
Metropolis updates have been replaced by a demon-update [29]. Most of the parameters 
in the algorithm were chosen by trial based on previous experience with the CP'^-model 
[25]. 

Our basic update steps are performed on 0(2) subgroups of the SO{N) group. We 
have chosen the same subgroup for each of the sites of the lattice. After a number of 
sweeps a new subgroup is chosen. 

The 0(2) subgroups that we consider are given by rotations among two rows or two 
columns. 

U'k^{x,t) = Si{x,t) Ukt{x,t) + S2{x,t) Ukj{x,t) 

U'kj{x,t) = -S2{x,t) Ukt{x,t) + si{x,t) Ukj{x,t) (3.1) 

with sl + si = 1. 

This parametrization induces an action for the embeded 0(2) model 

2 

^cond ('^) — ^ ^ ^ ^ ^mn {x ^t^ X , ^ ) Sjyi (x , ^) Sji (x,/^) (3.2) 

<x,t;x' ,t'> m,n=l 

with 

cii = C22 = J2^kz{x,t) Ukt{x',t') + Ukj{x,t) Ukj{x',t') 
k 

C12 = -C21 = J2^kz{x,t) Ukj{x',t') - Ukj{x,t) Ukt{x',t') (3.3) 

k 

For the updating of the embeded 0(2) model we used a combination of standard 
metropolis, demon-updates and microcanonical updates. We apply the microcanonical 
update step discussed by [27] for the standard XY model in two dimensions. First we 
compute the sum of the nearest neighbour spins of the site {x,t). 

2 

Rm= J2 ^ Cm„(a;,i,a;', i') s„(a;',i') (3.4) 

(x' ,t')nn(x,t) n=l 

The new values for s are then obtained by reflection with respect to R. 

s' = 2 5^R-s (3.5) 

Since Rs' = ils this update step keeps the action constant. 

The aim of the demon-update is to perform updates similar to a Metropolis update 
but avoiding CPU-intensive parts like the evaluation of trigonometric functions, the 
exponential function and pseudo-random numbers. 

The demons are introduced by an additional term in the action 

S' = S + J2dx,t (3.6) 
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applied to the demons. 

rf= -111(77), (3.7) 

where 77 is a pseudo-random number with a uniform distribution in the interval ]0, 1]. 
Then we perform updates that keep the composite action of the spin model plus the 
demons constant, exchanging energy between the demons and the spins. First we 
compute a proposal for a new spin- value s' by reflecting s at the sum of the upper and 
left neighbour spins. Then we check whether the demon at the site can take over the 
energy without becoming negative. If this is the case, we accept the proposal s' and 
set the demon to it's new value. After a sweep of such demon updates we translate the 
demons on the lattice. 

For one given embedding we performed 1 standard Metropolis sweep, and 10 to 
60 overrelaxation sweeps, and flnally 5 demon-updates. The number of overrelaxation 
sweeps was chosen to be roughly proportional to the correlation length. 

We alternate the row embedding and the column embedding. The flrst raw (column) 
is chosen in a flxed sequence from 1 to N, while the second was chosen randomly from 
the remaining ones. 

3.1 The improved correlation function estimator 

In order to obtain a meaningful result for the step-scaling function the correlation length 
on the flnite lattices has to be computed with an accuracy of less then one percent. In 
order to achieve this aim we used the improved estimator for the correlation function 
discussed in ref. [30] for the case of 0{N) vector models. 

The underlying physical idea for this improved estimator is similar as that for the 1- 
loop solution of the model discussed in section 2. However here instead of an effective 
one dimensional model we rather use a conditional (or embedded) model. For a given 
held conflguration U the conditional model is deflned by 

ScondiW) = -/3j2J2^r (W{t) U{x, t) {W{t + 1) U{x, t + 1))^) (3.1.1) 

t X 

where W is the held of the conditional model. Note that there is no spatial part in the 
action, since W{t) does not depend on x. Performing the x summation we get 

ScondiW) = -Y,Tr {Q{t,t+l)W^{t + l)W{t)) (3.1.2) 

t 

where 

Q{t,t+l) = -fiY, U{x, t) U^{x, t+1) . (3.1.3) 

X 

Reparametrising the model by R{t^t-\- 1) = W^{t)W{t-\- 1) we obtain 

Scond{R) = -J2Tr{Q{t,t+l) RT{t,t + l)) . (3.1.4) 

t 

For free boundary conditions in time-direction there is no constraint on the R^{t,t + 
1). Therefore the partition function factorizes, and the solution of the conditional 
1-D model is reduced to the solution of zero-dimensional systems. The conditional 
expectation value of the time-slice correlation function is given by 

{G{t, t + T))cond = Tr {S{t) {R{t, t + l))cond-{R{t + r - 1, i + r)),„„d S{t + r)) 

(3.1.5) 
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We were not able to compute the conditional expectation values exactly. Instead we 
used Monte Carlo integration for this task. Firstly one has to note that finite statistics 
for the conditional expectation value does not corrupt the end-result for {G{t,t-\- r)). 
Secondly a enormous gain in statistical accuracy can be obtained, since as a consequence 
of the factorization in eq. (3.1.5) the statistics of the single "baby" Monte Carlo's 
multiply. Typically we performed 200 Metropolis update steps for the evaluation of 
{R{t,t + l))cond- We used the final value of R for updating the field U. However, we 
have yet to make careful tests to determine the efficiency gain of this measurement 
technique. 



4 Simulation results for the 5*0(4) model 

The main object of this investigation was to test the predictions for ni/Aj^ for the 
SO{N) principal chiral models given by Hollowood [1, 20]. The action for the SO{N) 
matrix models is given in equation (1.1). It was found that even for moderate values of 
N {N > 6) the continuum limit was difficult to achieve with any degree of confidence, 
all indications being that a large correlation length is necessary. This is in contrast to 
recent work on SU{N) models [8, 9] where the results indicate that the models are close 
to the continuum limit even for small ^ (^ > 5). In the next section we will present 
some results for = 6, 8 and 10 to support these statements but will postpone our 
speculations concerning why such a difficulty occurs until the conclusion. 

In this section we concentrate on the 5*0(4) model and apply the renormalization 
scheme described in [2], the '.^—scheme', and compare it with the lattice- and energy- 
schemes used in [8, 9]. The coupling constant u{Lmoo) is defined in equation (1.3) 
[24] and to achieve the continuum limit we require u{Lmoo, a/ L) for fixed Lnioo to be 
essentially independent of a and hence that ^oo is large enough that the a-dependent 
corrections to equation (1.3) are negligible. The problem is that in practice this might 
require ^oo to be very large indeed and unachievable in a present-day simulation. Alter- 
natively, it might be possible to fit the a— dependence in equation (1.3) by measuring 
u{Lmoo, a/ L) for various ^oo with Lnioo fixed, and extrapolating to a/L = . In the 
case that the corrections are perturbative they behave as 0{a^) . Whether or not this 
is the case must be deduced from the simulation and for the 5*0(4) model there is clear 
evidence that a simple perturbative interpretation of the a— dependent effects is not 
possible for the values of ^oo we use. Nevertheless we try to extrapolate the results to 
the continuum limit in a reliable way and compute the value of u{2,a/L). It must be 
emphasized that it is crucial to determine this value of u as accurately as possible since 
it is the starting point for the subsequent determination of u at smaller scales, and 
ultimately contributes to the systematic error on the computation of ni/Aj^. In order 
to attain the continuum limit we measured u{2) for increasing values of ^oo and fitted 
the a— dependence and so deduced the continuum coupling, u{2). We characterised the 
theory either by the lattice coupling (3l or by the coupling defined by Parisi [19] in 
terms of the internal energy, equation (2.1.1): 

- 1 . ^ 

Note (E) is normalized to lie in [0, 1] . We found that (E) suffered from strong finite-size 
effects and so could not be expressed as a function of (3l alone. We chose to measure 
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by trial we can see which method gives the better extrapolation to a/L = 0. Where 
L/a = 2^00 was non-integral, simulations were done for the integer L/a values either 
side of 2^00 and interpolation used to deduce S,{L) ■ The results are shown in table 1 
where we can see that the energy-method shows better convergence than the lattice- 
method. The errors shown are statistical and a naive straightforward extrapolation 
gives the value u{2) = 2.25(2) for the continuum coupling constant. However, 
we shall argue below that whilst plausible this value for u{2) is incorrect since the 
convergence to the continuum is only apparent not real. 



h 


(E) 


Coo 


u(2,a/L) 


u(2,a/L) 








lattice-method 


energy-method 


1.05 


0.5411(2) 


3.71(3) 


1.98(1) 


2.44(1) 


1.10 


0.4702(2) 


8.34(7) 


2.07(1) 


2.35(2) 


1.12 


0.4435(1) 


13.62(13) 


2.11(1) 


2.27(2) 


1.14 


0.42176(6) 


25.3(4) 


2.17(2) 


2.26(2) 



Table 1: The values of u{2,a/L) in lattice- and energy-methods as a function of ^00 • The 
energy-method shows plausible convergence to the extrapolated value u{2) = 2.25(2) . 



That (E) is sensitive to the lattice size is another indication that SO{N) models 
are more complex in the continuum limit than SU{N) models where no such effect 
is observed. We find that the finite-size effects in (E) largely seem to offset those in 
u{Lmoo) which indicates that there is some connection between the short and long- 
range properties of the system. From table 1 it can also be seen that the a— dependent 
effects in the lattice-method are not fitted by a perturbative parametrization: they are 
closer to a a^/^ dependence. This indicates that non-perturbative contributions are 
strong and that their effect on the calculation of m(2) seems to be largely accounted 
for using the energy-method. These are all reasons why we should be cautious and 
suspicious of assuming that we are observing properties of the continuum theory. 

One reasonable test that we are simulating the continuum theory is to check that the 
computed value of m/ hjjg agrees with the theoretical prediction. We thus compute the 
step-scaling function described in [2] and section 1 and fit the short-distance behaviour 
of u to the form deduced from perturbation theory. We choose the factor for the 
scale-change to be s = 1/2. Therefore we consider pairs of lattices with sizes 2L and L 
respectively and /3/, is adjusted so that u{2L/a, (3l) is a required value. Then u{L/a, (3l) 
is measured in both the lattice and energy schemes, that is, keeping either (3l or (3e 
constant on the lattice pair. Various values of L/a were chosen so that the a— dependent 
effects can be determined and eliminated by extrapolation. When we believe that the 
continuum limit has been attained then the step scaling function a can be determined 
from equation (1.5) and compared with perturbation theory using equations (1.6,1.7) 

^° - 27r(7V-l) - y • 

The first two coefficients of the /3— function quoted for u are obtained from those asso- 
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lattice-method 


energy-method 


8 


1.0908 


0.4531(5) 


2.25 


1.590(10) 


1.727(10) 


14 


1.1169 


0.4352(3) 


2.25 


1.647(8) 


1.720(8) 


26 


1.1393 


0.4186(2) 


2.25 


1.711(7) 


1.737(7) 


40 


1.15335 


0.40915(7) 


2.25 


1.734(7) 


1.747(7) 


7 


1.1088 


0.4347(3) 


1.747 


1.416(3) 


1.475(3) 


13 


1.1347 


0.4194(2) 


1.747 


1.469(3) 


1.480(3) 


20 


1.15235 


0.4082(1) 


1.747 


1.492(3) 


1.499(3) 


30 


1.16917 


0.39778(6) 


1.747 


1.506(3) 


1.509(3) 


10 


1.1484 


0.4111(1) 


1.517 


1.323(2) 


- 


15 


1.1659 


0.3989(1) 


1.517 


1.340(3) 


- 


20 


1.1803 


0.3911(1) 


1.517 


1.345(2) 


- 


5 


1.1424 


0.4074(3) 


1.351 


1.203(2) 


- 


10 


1.1789 


0.3906(2) 


1.351 


1.212(2) 


- 


16 


1.2045 


0.3784(1) 


1.351 


1.218(2) 


- 


8 


1.2028 


0.3774(2) 


1.222 


1.113(1) 


- 


16 


1.245 


0.36014(9) 


1.222 


1.114(1) 


- 


8 


1.2452 


0.3587(2) 


1.115 


1.025(1) 


- 


16 


1.2914 


0.34215(8) 


1.115 


1.024(1) 


- 


8 


1.2922 


0.3408(2) 


1.0234 


0.9494(7) 


- 


16 


1.3368 


0.32681(7) 


1.0234 


0.954(1) 


- 


8 


1.3357 


0.3264(1) 


0.948 


0.8902(5) 


- 


16 


1.3851 


0.31230(6) 


0.948 


0.8910(4) 


- 


8 


1.3849 


0.3116(1) 


0.891 


0.8355(5) 


- 


16 


1.4356 


0.29861(6) 


0.891 


0.8342(5) 


- 



Table 2: Couplings measured in the lattice- and energy-methods for different values of L/a. 
The a-dependent violations to finite-size scaling are most pronounced for Lnioo ~ 1. The 
energy-method was not used where it gave little information in addition to the lattice-method. 
The statistical errors oi u[2L / a) which are not quoted here are of similar size as those of u{L/ a). 
coupling values. 



elated with ul = 1//3l [21] by using equation (2.1.25) gives the tree-level relation 

{N-l) 

The results are shown in table 2. 

In each case the continuum coupling for a given scale was determined by extrapola- 
tion to a = . From table 2 it can be seen that the most care needs to be taken when 
Lnioo ~ 1 where the corrections to finite-size scaling are the greatest. In all cases the 
energy-method was the most convergent, but for sufficiently large L both the lattice- 
and energy-methods were compatible and a common value for the continuum coupling 
was consistent. For Lnioo < 1/2 the a— dependent corrections were small and consistent 
with perturbation theory and both schemes gave consistent results. Following eq. (1.8) 
we tried to use the result of a given step as the argument of the next step. This was 
achieved in all cases except when Lnioo = 1/2 and Lnioo = 1/128 where twe corrected 
the small mismatches by interpolation. These corrections lead to scale changes of 2.026 
and 1.840 respectively. The sequence of continuum couplings deduced are shown in 
table 3. 
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Figure 1: Values of the coupling u{Lmoo,a/L) from table 2 plotted against L/a Lnioo fixed 
at the values 2, 1, 1/2 . The violations of the finite-size scaling assumption are clearly evident 
for both versions of the ^—scheme used: the lattice-scheme (•) and the energy-scheme (■). In 
each case the value assumed for the extrapolation to a = is shown as the dashed line 

Also shown are the 1— and 2— loop approximations to a{l/2,u) evaluated with 
u = u{2Lmoc) ■ These should be compared with the Monte-Carlo result 



Clearly asymptotic scaling sets in for Lnioo < 1/4 . It is perhaps surprising that 
the 2— loop approximation fits the Monte-Carlo results so well at such relatively large 
scales, but it gives confidence that we can probe deeply into the region where asymptotic 
scaling holds and hence that perturbative parametrization of u in terms of is valid. 
The a-dependence of the results for these smaller values of u indicates that (3l is large 
enough for the extrapolation to the continuum limit, a = 0, to be reliable. The data 
from table 3 are shown in figure 2 where the computed function a{2, u) is compared with 
its 1— and 2— loop approximations deduced from equations (1.6,1.7). The important 
intermediate region where it is crucial to maintain the continuum limit in order to 
reliably relate the low and high energy scales is also in clear evidence. This match 
between the small-L region and large-L region (where trivially u{Lmoo) = Lrrioo) is the 
important part of the simulation. 

From table 3 we can deduce the value of ni/A^, i.e., in the ^—scheme, using the 
2— loop formula, equation (1.9). We have: 



(t(1/2, M(2Lmoo))MC = u{Lmco) ■ 



1 




(4.2) 
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i-loop 


/-loop 


rf-loop 


1 


2.25 


1.754(13) 


1.93054 


1.9010 


1.8915 


.987/2 


1.747 


1.517(6) 


1.5481 


1.5329 


1.5291 


.987/4 


1.517 


1.351(6) 


1.3647 


1.3544 


1.3521 


.987/8 


1.351 


1.222(4) 


1.2289 


1.2214 


1.2199 


.987/16 


1.222 


1.1148(21) 


1.1212 


1.1154 


1.1145 


.987/32 


1.1148 


1.0234(15) 


1.0303 


1.0259 


1.0252 


.987/64 


1.0234 


0.9478(11) 


0.9518 


0.9483 


0.9478 


1.073/128 


0.9561 


0.8912(12) 


0.8933 


0.8904 


0.8900 


1.073/256 


0.8912 


0.8339(12) 


0.8364 


0.8340 


0.8337 



Table 3: The sequence of continuum couplings as a function of Lnioo ■ Also shown are the 
1— and 2— loop approximations to a{l/2,u) evaluated with u = u[2Lmco). These should be 
compared with the Monte-Carlo result (t(1/2, M(2_Lmoo))MC = u[Lmco). It seems plausible 
that asymptotic scaling sets in for Lnioo < 1/4 . 



In fact, since the covering group of 5*0(4) is SU{2) ® SU{2) and the SU{2) matrix 
model is isomorphic to the 0(4) spin model, we are able to use the 3-loop /3-function 
in reference [2] to deduce that 

A discussion of the detailed relationship between these models will be postponed until 
the next section but it is convenient to invoke this formula here. The 1-loop calculation 
necessary to determine the ratio A^/Aj^ is given in section 2, and using the result 
for Al/A^ and Aj^/Al from equations (2.1.4) and (2.1.23), we give the computed 2- 
and 3-loop results for ni/Aj^ versus length scale in table 4. The prediction for this 
quantity can be taken from the paper by Hollowood [1] where the formula for ni/Aj^ 
for the SO{N) matrix models can be extended down to = 4 [20]. The prediction is 

m 2^/2 ^ ^ 

: 3.8715 . (4.4) 



Lnioo 


1/256 


1/128 


1/64 


1/32 


1/16 


1/8 


1/4 


1/2 


1/1 


2-loop 


14.3(1.0) 


14.3(1.0) 


14.4(9) 


14.4(9) 


14.1(9) 


14.1(9) 


14.2(9) 


14.0(7) 


13.2(5) 


m/Aj^ 
3-loop 


13.7(1.0) 


13.6(1.0) 


13.7(9) 


13.6(9) 


13.3(9) 


13.2(9) 


13.2(9) 


12.9(7) 


12.0(5) 



Table 4: Computed values for m/ A-j^ using the 2-loop /3— function for the coupling evaluated 
at the indicated scale (the small adjustments to the scale explicit in table 3 are omitted for 
clarity). The errors are calculated from the accumulation of statistical errors at all preceeding 
steps. 



As can be seen from table 4 the computer results are wrong by a minimum of a 
factor of four. The inescapable conclusion is that we have not eliminated corrections to 
finite-size scaling and are not close to the continuum limit in our simulation. Without 
the theoretical prediction we might have been persuaded that the evidence did point to 
the simulation having produced reliable continuum results - the step scaling function 
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strongly violated is for lattices where Lnioo ~ 1 and the parts of the simulation which 
are important for setting the scale are not probing the continuum. The corollary is that 
we need the inital value of ^oo = ^/arrioo to be considerably larger than the biggest we 
have taken. The surmise is that since the SO{N) manifold is not simply connected, 
Ili{SO{N)) = Z2, there are vortices which are responsible for non-perturbative viola- 
tions of finite-size scaling and much larger values of (3 are needed before their effect is 
sufficiently suppressed so that the continuum limit can be approached in a controlled, 
perturbative, way. One check is therefore to simulate the covering group of 5*0(4), 
SU{2) ® SU{2), and see whether the problems with finite-size scaling violations are 
perturbative. The results are presented in the next section. 

We can also check what happens if we apply the method used in references [12, 10, 
8, 9] for the couplings defined in equation (1.11). These bare couplings are defined to 
be evaluated at the scale of the lattice spacing a. The A-parameters of the respective 
lattice- and energy-schemes, Al and can be calculated using equation (1.9). Since 
the corresponding value of anioo is known from the simulation the estimate for ni/Aj^ 
can be calculated using equations (2.1.3) and (2.1.4). The results are shown in table 5. 
The results are very poor and it is clear that the approach is hopeless. This should be 
contrasted with the success of the method applied to SU{N) matrix models [8, 9] and 
the 0(4) and 0(8) spin models [10]. 



h 


(E) 


Coo 




lattice-scheme 


energy-scheme 


1.05 


0.5411(2) 


3.71(3) 


305(2) 


20.4(2) 


1.10 


0.4702(2) 


8.34(7) 


248(2) 


31.3(3) 


1.12 


0.4435(1) 


13.62(13) 


194(2) 


34.1(3) 


1.14 


0.42176(6) 


25.3(4) 


133(2) 


30.9(5) 


1.15 


0.41291(5) 


33.4(8) 


114(3) 


29.5(7) 



Table 5: The values of m/ A-j^ in the lattice- and energy-schemes [12, 10, 8, 9] as a function 
of 'Coo for the 50(4) model. The theoretical prediction is m/A-j^ = 3.87153 [1, 20] 



5 Results for the SU{2) (g) SU{2) covering of 5*0 (4) 

The SU{2) ® SU{2) matrix model consists of two independent SU{2) models and the 
SU{2) matrix model is isomorphic to the 0(4) spin model where the fields take values 
in 0(4)/0(3) = S'^. Thus we need only to simulate the 0(4) spin model which can be 
done using the cluster algorithm [11]. The theoretical prediction for m/ hjjg is given 
in [7, 6] to be 

WAm5 = = ^-^^^^ • (^-^^ 

The only difference between the 5*0(4) and SU{2) matrix models is that the lightest 
state of the SU{2) model is a spinor and that of the 5*0(4) model corresponds to a 
spinor-antispinor state which is not, in fact, bound. Thus, in the 50(4) model the large 
time asymptotics of the correlator are controlled by the spinor-antispinor cut and not 
a bound state pole. However, in two dimensions the large time behaviour is dominated 
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1^30(4) = 2m5C/(2) 



(5.2) 



This factor of two simply converts the 5*0(4) prediction of Hollowood [1] into the 
prediction of Hasenfratz et al. [7, 6] for the 0(4) spin model. It then follows that the 
two continuum couplings defined as in equation (1.3) are related by a factor of two: 

^^50(4) = 2m5[/(2) . (5.3) 

We shall omit the distinguishing subscript on u unless it is necessary to avoid ambiguity. 
The action for the 0(4) spin model is taken to be 

^(s) = /3l^s„-s„+^ . (5.4) 

As before we define u using equations (1.2) and (1.3). The tree-level result from 
section 2.2 relating the ^- and lattice-schemes is 

= . (5.5, 

Using this result, the /3— function for u as defined in equation (1.7) then has coefficients 
. _ (A^-2) ^ _ N-2 ^ N-2 



7r(7V- 1) ' ' 7r2(7V- 2)2 ' ' 7r^{N 



(5.6) 



The result for 62 is taken from reference [2]. In addition to the 2-loop formula for A^, 
equation (1.9), we have the 3-loop formula 

,(3) - a(2) , , _ 

k{N - 1) 



Af) = A'f)[i-—^u) . (5.7) 



Because the continuum limits of the 5*0(4) and 0(4) spin model are controlled by the 
same Lie algebra the conversion ratio A^/Aj^ is the same for both. However, as a 
check and for completeness, the 1-loop calculation which yields this conversion ratio 
for the general O(A^) spin model is briefiy described in section 2.2. Of course, for = 4 
all necessary results can be taken from reference ([2]). 





(E) 


X 


Coo 


«(1, 1/^00) 


2.00 


0.4230(2) 


99.8(3) 


7.89(2) 


1.552(5) 


2.20 


0.3775(1) 


266.6(8) 


13.95(5) 


1.562(4) 


2.40 


0. .34070(5) 


555.7(1.9) 


25.7(1) 


1.585(5) 


2.60 


0.31071(2) 


1629.9(2.6) 


47.07(8) 


1.584(3) 



Table 6: The values of u{l, l/^oo) in the lattice-method as a function of ^00 . Interpolation is 
used to calculate u at non-integer values of L. Clearly, finite-size scaling violations are small 
and, within the errors, not incompatible with the perturbative prediction that they behave like 

The simulation results for u{l) for various values of ^00 are shown in tables 6 and 7 
and the results for the corresponding step-scaling function are given in table 8 where 
the values quoted have been extrapolated to the a = limit. Clearly the violations of 
finite-size scaling are much smaller than for the 5*0 (A) matrix model, compare with 
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7 


2.1914 


0.3714(2) 


47.87(6) 


1.584(4) 


1.210(3) 


13 


2.4032 


0.3378(1) 


140.8(2) 


1.584(3) 


1.219(2) 


24 


2.6050 


0.30941(3) 


415.3(4) 


1.584(3) 


1.227(1) 


8 


2.4649 


0.3264(1) 


78.0(1) 


1.288(3) 


1.008(2) 


16 


2.6959 


0.29735(5) 


268.0(3) 


1.288(2) 


1.011(2) 


8 


2.6954 


0.29579(7) 


96.5(1) 


1.011(2) 


0.865(1) 


16 


2.9260 


0.27185(4) 


337.9(4) 


1.011(2) 


0.863(1) 



Table 7: Couplings measured in the lattice-method for different values of L/a. The a-dependent 
violations of finite-size scaling are apparent but small, even for the crucial case where L/a ^ ^oo 



Lnioo 


u{2Lmoo) 


u{Lmoo) 


a(l/2,u(2Lm^)) 


1-loop 


2-loop 


3-loop 


2/1 


4.132(10) 


2.309(10) 


2.5700 


2.2927 


2.1062 


1/1 


2.309(10) 


1.584(4) 


1.7236 


1.6384 


1.5925 


1/2 


1.584(4) 


1.228(3) 


1.2847 


1.2495 


1.2348 


1/4 


1.228(2) 


1.011(3) 


1.0401 


1.0216 


1.0152 


1/8 


1.011(2) 


0.863(2) 


0.8801 


0.8689 


0.8657 



Table 8: The sequence of continuum couplings for the 0(4) spin model as a function of 
Lnioo ■ Also shown are the 1— and 2— loop approximations to a{l/2,u) evaluated with u = 
u{2Lmoo)- These should be compared with the Monte-Carlo result a{l/2, u{2Lmoo))MC = 
u{Lmoo)- Clearly, asymptotic scaling is already setting in for Lnioo < 1 

tables 1 and 2, and are compatible, within errors, with the perturbative prediction 
that they behave like 0{a^). In comparing tables 7 and 2 the conversion factor of two 
between couplings, eqn (5.3) should be bourne in mind. Note that we have not used 
the energy-method for the spin model since it is not needed. 

The 2-loop and 3-loop computed values of ni/Aj^ are given in table 9. The 2- 
loop result is already near to the predicted value of 1.9358 for Lnioo ~ 1/8 and the 
3-loop result agrees with this prediction within errors even for Lnioo ~ 1. In figure 2 
we compare the SU{2) ® SU{2) and 5*0(4) results for the step-scaling function. We 
have plotted uso(4) and 2usu(2) on th^ ordinate since in the continuum limit the data 
points should coincide up to the small correction factors of table 3. Also shown are the 
curves for a{l/2, u) derived from the 1-loop and 3-loop /3-functions. The large u result, 
(t(1/2, m) —7- u/2, to which all curves should eventually be asymptotic is also shown. It 
can be seen that there is a clear deviation between the results for the two models in the 





1/8 


1/4 


1/2 


1/1 


m/Aj^ 
2-loop 


1.783(30) 


1.735(28) 


1.679(18) 


1.609(10) 


3-loop 


1.963(33) 


1.944(31) 


1.930(20) 


1.935(11) 



Table 9: Computed values for m/ Kjjg using the 2-loop and 3-loop /3— functions for the coupling 
evaluated at the indicated scale. The errors are calculated from the accumulation of statistical 
errors at all preceeding steps. These results are to be compared with the theoretical prediction 
[7, 6] m/A^= 1.9358 
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Figure 2: The step-scaling function (t(1/2,m) versus u. The solid and dashed curves are 
respectively the 1- and 3-loop perturbative calculations. Data for the 50(4) matrix model (o) 
and the SU (2) ®SU (2) matrix model (□). The large u result to which all curves are asymptotic 
is shown as the long-dashed line. 

The calculations of m/Aj^'m the lattice- and energy-schemes have been done using 
3-loop results by Wolff in ref. [10]. Our simulation results agree in detail with his and 
he finds general agreement with the prediction for ni/Aj^. However, the different 
schemes tried by Wolff do show different rates of convergence to theory as a function of 
^oo the best being the energy-scheme which agrees very well with theory for ^oo ~ 10. 
We have similarly good agreement for the ^-scheme confirming the ease with which the 
continuum limit can be controlled. 

6 Results for SO{N), iV = 6, 8, 10 

We attempted to analyse the SO{N) models for = 6, 8, 10 in the same way as for 
5*0(4). However, the CPU time required is prohibitively large and we were unable to 
work with sufficiently large correlation lengths. It is instructive, however, to compute 
the values of m/Aj^for each model in the lattice- and energy-schemes. The results are 
shown in table 10. Whilst the energy-scheme does not show the large deviation from 
theory of the 5*0(4) model the results are clearly untrustworthy. There is no agreement 
between the different schemes and, although the range of ^ used was limited, there is 



20 



for = 8, 10 the energy-scheme gives results which seem independent of ^oo for the 
restricted range covered. Clearly, this cannot be taken as indicating that the results 
have converged to the ^oo co limit: in light of our experience it shows very little. 

Another prediction derived from the exact S'-matrix is for the mass-ratios of particles 
in the theory. For SO{N) {N even) the prediction is 

m„ sin(7r»/(A^ — 2)) , ^ , , ^ 

—^ = —\J^ — -11 l<p<7V-2/2, 6.1 
m sin(7r/(7V - 2)) ^ " ^ " ' ^ ' 

where p labels the p-th species and m is the mass of the lightest state. We simulated 
the SO{N) model for = 6, 8, 10 and measured the masses in the different channels 
labelled by p. Hollowood [1] has discussed the relevant interpolating operators for these 
states and we choose the simplest operators which couple to the desired state in each 
channel. The operator for the p-th state is 

0,(x).,...,,,...,, = Yl ' (6-2) 

perms of ji...jp 

where P is the permutation signature of the ordering of the {ji}. Thus Op is the outer 
product of p matrices antisymmetrized on the row and column labels respectively. The 
corresponding Green function is 

^.(x) = (rr(O,(x)Oj(0))), , (6.3) 

where the trace has the obvious meaning. The results for the mass ratios 1712/171 for 
A^ = 8, 10 and m-^jm for A^ = 10 are shown in table 11. 



N 


h 


(^) 


•Coo 




lattice-scheme 


energy-scheme 


theory 


6 


1.5 


0.6786(1) 


1.387(3) 


48.8(1) 


5.32(1) 


3.87153 


6 


1.6 


0.6226(1) 


1.97(1) 


62.3(3) 


6.03(3) 


3.87153 


6 


1.7 


0.5181(3) 


4.85(4) 


46.0(4) 


7.98(7) 


3.87153 


6 


1.73 


0.4784(2) 


8.7(1) 


30.7(3) 


8.07(9) 


3.87153 


8 


2.28 


0.5727(3) 


2.80(2) 


35.3(3) 


5.27(4) 


3.65837 


8 


2.31 


0.5362(4) 


3.98(3) 


28.0(2) 


5.55(4) 


3.65837 


8 


2.33 


0.5073(3) 


5.75(7) 


21.0(3) 


5.51(7) 


3.65837 


8 


2.35 


0.4851(2) 


8.05(6) 


16.2(1) 


5.36(4) 


3.65837 


10 


2.91 


0.5623(2) 


3.05(1) 


25.13(8) 


4.65(2) 


3.523789 


10 


2.93 


0.5341(4) 


4.18(3) 


19.46(14) 


4.61(3) 


3.523789 


10 


2.95 


0.5054(2) 


5.90(3) 


14.63(7) 


4.62(2) 


3.523789 



Table 10: The values of mj h-^ computed in the lattice- and energy-schemes compared with 
the theoretical prediction for N = 6,8, 10. There is no agreement between the schemes and no 
trend suggesting that the results will converge to the prediction. 

There is no convincing agreement between simulation and theory and, moreover, 
no trend suggesting that the discrepancy is ~ a^. 



7 Discussion 



The main result of this paper is that the properties of the continuum 5*0(4) theory 
cannot be observed in a simulation of the lattice-regularized model for the values of /3/, 
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simulation 


tiieory 


simulation 


tneory 


8 


2.28 


2.80(2) 


1.84(1) 


1.732 


- 


- 


8 


2.31 


3.98(3) 


1.79(3) 


1.732 


- 


- 


8 


2.33 


5.75(7) 


1.83(2) 


1.732 


- 


- 


8 


2.35 


8.05(6) 


1.85(2) 


1.732 


- 


- 


10 


2.91 


3.05(1) 


1.89(2) 


1.848 


2.7(2) 


2.414 


10 


2.93 


4.18(3) 


1.93(2) 


1.848 


2.73(3) 


2.414 


10 


2.95 


5.90(3) 


1.91(2) 


1.848 


2.73(3) 


2.414 



Table 11: The computed mass-ratios for N = 8, 10 compared with the predictions. There is 
no convincing agreement between simulation and theory. 

and correlation lengths accessible to current computers. We have shown that no such 
difficulty occurs for the model based on the SU{2) ® S'f/(2) cover of 5*0(4). We believe 
that both models give rise to the same continuum theory, characterized by the fixed 
point at (3l = oo, but that the ways in which this continuum theory is approached in the 
lattice-regularized versions are very different. For SU{2)(SiSU{2) finite-size scaling with 
0{a^) deviations holds for the range of couplings used, and for scales Lnioo < 1 the flow 
with L of the renormalized coupling, u{Lmoo), is well given by the 3-loop /3— function. 
The value for m/ hjjg computed in the simulation agrees well with the theoretical 
prediction (table 9). In contrast, for 5*0(4) the violations of scaling do not fit an O(a^) 
form but do seem to diminish to zero (see fig. 1) as ^oo increases. This apparent or 
'pseudo-scaling' can be wrongly interpreted as signalling the continuum theory and can 
deceive us into believing that the value for the continuum coupling can be deduced. 
The false nature of this pseudo-scaling is exposed by comparing the resulting computed 
value for m/ of ~ 14 with the theoretical prediction of 3.8715. We also find no 
convergence to a result for m/hjjg using the lattice- and energy- schemes for 5*0(4). 
In contrast for SU{2) ® SU{2) the lattice-scheme gives an acceptable result although 
clearly inferior to that of the ^—scheme. 

We analysed data for the SO{N) models with = 6, 8, 10 and used the lattice- and 
energy-schemes to attempt to obtain an estimate for mhjjg. The results are shown in 
table 10 and it is clear that there is no agreement between the schemes nor with the 
theoretical prediction. We also computed the mass-ratios of the fundamental masses 
predicted by Ogievetsky et al. [5] and the results are given in table 11. There is a 
persistent discrepancy up to 10% and there is no sign of a trend to the correct values 
as [iL increases. These results support the conjecture that we are unable to simulate 
the continuum theory for SO{N) models with present computer resources. 

The discrepancy between the computed and theoretical values of m/ hjjg is much 
larger than that of about 10% reported by Liischer et al. in their analysis of the 0(3) 
spin model [2], which they attributed to the truncation of the perturbative /3— function 
at 3-loop order. This is not the resolution of the problem we have found for the 
50(4) model. We conjecture that the difference lies in the different connectivities of 
the underlying manifolds: the SO{N) models contain Z2 vortex lattice artifacts whilst 
the covering group models do not. In the general SO{N) case the cover is Spin{N) 
which is constructed from the associated Clifford algebra [31]. The vortices create an 
obstruction to observing the fixed point at /3/, = 00 but are eventually suppressed at 
sufficiently large (3l. Recent work by Hasenbusch [16] and Niedermayer et al. [17] 
has discussed a similar phenomenon comparing the 0(3) and RP^ spin models. They 
propose a similar conclusion, namely that the difference between the two models is due 
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deduce that u{2) = 2.25 but from our simulation of the covering group we find that 
the value should be m = 3.17. We expect that as (3l increases a cross-over phenomenon 
will occur where the violations to scaling will again become large and then eventually 
diminish to become 0{a^) allowing the true scaling limit to appear. It is important 
to estimate the value of (3l and ^oo above which finite-size scaling and the continuum 
theory limit should be observed. A crude attempt can be made with current data by 
using the theoretical prediction for ni/Aj^ together with 

= -^{bouiLm^yy/'eMj-^ ,) , (7-1) 

m/A^ bou{Lmoo) 

to deduce ^oo{Pl) from the data for u{Lmoo) at sufficiently small Lrrioo given in table 3. 
These results for ^oo can be compared with the results for ^oo computed directly 
from simulation. If the lattice theory is near the continuum limit these alternative 
methods of computing ^oo{Pl) should give similar answers. There is consistency to 
within 10% for the values ^oo{Pl) from using similar values of (3l on lattices of different 
widths corresponding to different measured values of u. Where there is a choice we 
have taken the result from the smallest u. We plot log(^oo) versus (3l in figure 3 for 
both approaches. The mismatch is clear for (3l ~ 1.15: the direct measurement gives 

= 33.4(8) whereas the short-distance result is ^oo ~ 150. It seems reasonable to 
infer that the two methods will not agree until ^oo ^ 150. 

The perturbative step-scaling function is well reproduced for sufficiently small scales 
even though the observables on scales Lnioo ~ 1 show a large departure from continuum 
behaviour. In this context note that a{l/2,u) for u = 2.25 is actually quite close to 
the 3-loop prediction even though the true scale associated with this value is a factor 
of about four different from that assigned in the simulation. Thus near agreement with 
the perturbative prediction at small scales is not sufficient to infer that the continuum 
theory is being observed at large scales: from our simulation of the covering group the 
correct value for u at the scale assigned in this case is m = 3.17 and a{l/2, u) for this 
value agrees very well indeed with the 3-loop perturbative prediction (see fig 2). 

Suppose we were able to simulate at, say, (3l ~ 1.18 on a large enough lattice. From 
our results we see that the properties of the continuum theory are well reproduced at 
scales Lnioo < 1/8 but from the above discussion we would expect the direct measure- 
ment of ^oo to be considerably less than the short-distance prediction of ^oo ~ 300. So 
whilst the properties of the continuum theory are computable at short distances for 
(3l ~ 1-18 the long distance results do not refiect the continuum but are dominated by 
residual lattice artifacts: vortices in the 5*0(4) model. 

For an euclidian continuum field theory a field configuration can be viewed as a 
map of the space onto the manifold of the field. A vortex is now characterized 
by the property that the map of a loop in onto the manifold of the field is not 
smoothly contractible. As a consequence there is at least one sigularity of the field 
inside such a loop. In statisitical mechanics, vortices have been mainly discussed in 
relation with the two dimensional XY model. The classical energy of a vortex is given 
by ~ TT log (i?/a) where R is the size of the vortex. Based on the simple energy versus 
entrophy argument that the free energy is given hj F = E — TS with S = 21og(i?/a), 
Kosterlitz and Thouless [33] inferred the occurrence of a phase transition at T = 7r/2. 

This argument does rely on the assumptions, that the free energy of a vortex at a 
fixed location is essentially given by its energy and that knowledge of the free-energy 
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a fixed location is bounded as i? — ?■ oo for non-abelian tlieories [32] and so a simple KT 
style analysis cannot be carried out. However, it has also been suggested [34, 32] that 
for non-abelian theories the interaction between vortices is such that the free-energy of 
multi-vortex configurations cannot be inferred from the properties of an isolated vortex. 
It is clear that further studies are necessary to clarify the true position. 
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Figure 3: log(^oo) plotted versus /3i derived from the short-distance behaviour and the 2-loop 
/3-function (•) and from direct measurement in the simulation (■). The two sets of results do 
not agree indicating that the large scale properties are not controlled by the continuum theory. 

A number of approaches can be taken: 

(i) The free-energy of an isolated vortex can be calculated at least in 1-loop pertur- 
bation theory and in a lattice simulation to check that the argument above can 
be made. Also the free-energy of multi-vortex configurations should be computed 
by simulation; 

(ii) The vortices of SO{N) can be suppressed by eliminating any configuration which 
contains one or more vortices. The finite-size scaling analysis can be repeated to 
see whether the continuum theory is more readily observed; 

(iii) The matrix models based on Spin{N) can be simulated and compared with the 
SU{N) models and the mass-ratios predicted from the exact S'-matrix can be 
computed; 
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S{s^ s) = /3 ^ ai{x^t)s{x^t) • s{x + 1, ^) + a2{x^t)s{x^t) • s{x^t + 1) 

x,t 

- log(^) Pi'^) ' (7-2) 

plaquettes 

where s is an A^— component vector of unit length, ci^, = 1, 2, is a gauge 
field taking values in [1, —1] and P{cf) signifies the plaquette. 

(b) The model whose manifold has the same topology as the 5*0(4) manifold 
but in which the vortex operators can be explicitly constructed. This model 
has action 

S'((T,s,r) = /3 ^ (Ti (s, (s(a;, • s(a; + 1, + r(a;, • r(a; + 1, i)) 

x^t 

+ (J2{x^ t) (s(a;, t) ■ s{x, i + 1) + s(a;, t) ■ s{x, t + 1)) 
- log(^) Yl ^(^)' (7-3) 

plaquettes 

where s and r are four-component vectors of unit length and ci^ and P{(j) 
are as defined in (a). 

In both models the renormalization group fiow in (/3, z) can be studied using 
Monte-Carlo methods and hence the effect of vortices, measured by P{(j) 
and controlled by the fugacity, z, can be determined. 

These projects are currently in hand. 

The outcome is that we should be wary of claims that the continuum theory has 
been observed which are based on the observation of scaling in a limited window in the 
coupling constant. Even if properties of the continuum are observed at short distances 
it does not follow that observables on the scale of the correlation length, which are 
sensitive to so-called non-perturbative effects, are controlled by the continuum theory. 
This could be the case for any lattice model which has non-trivial topological artifacts. 
It has been pointed out [35] that QCD is such a theory since the gauge group is 
SU{3)/Z3 where is the centre of SU{3) [36]. It is important to determine whether 
such an effect exists in QCD and at what level of accuracy it needs to be taken into 
account in present-day simulations. We have been unable to observe any continuum 
properties in the SO{N) matrix models but it remains to be seen whether Spin{N) 
models have the same problems. 
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